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Abstract 

Weights of edges and nodes on food webs which are available from the empiri- 
cal data hide much information about energy flows and biomass distributions 
in ecosystem. We define a set of variables related to weights for each species 
i, including the throughflow Tj, the total biomass Xi, and the dissipated flow 
Di (output to the environment) to uncover the following common patterns in 
19 empirical weighted food webs: (1) DGBD distributions (Discrete version 
of a Generalized Beta Distribution), a kind of deformed Zipf's law, of energy 
flow and storage biomass; (2) The allometric scaling law Tj oc X®, which 
can be viewed as the counterpart of the Kleiber's 3/4 law at the population 
level; (3) The dissipation law Di oc Zf ; and (4) The gravity law, including 
univariate version oc (TjTj) 7 and bivariate approvement /y oc T^Tj 2 . 
These patterns are very common and significant in all collected webs, as a 
result, some remarkable regularities are hidden in weights. 

Keywords: Weighted food webs, Energy flow distribution, Scaling 
relations, Allometric scaling, Gravity law, DGBD rank-ordered curve 



1. Introduction 



Complex network is a useful tool to study interactions and re l ation- 
ships between components of a complex system ( j Watts and Strogatzl . Il998 ; 
Albert and Barabasil . 120021 ) . In ecology, complex network models a re used in 



many ways, including representing trophic interactions in food websdMontova and Sold. 



2002 



3D 

Dunne et all 120021 ; [Williams et all 2002; Berlo w et aUl2004HEmmerson and R affaelli. 



2004 ) and energy- matter flux in ecosystems(Odum . 1988 ; Finnl.ll976l;ISzvrmer and Ulanowicz . 



19871 ; lHigashil . ll986uHigashi et al.L Il993l ; iBaird and Ulanowiczl . ll989uFath and Patten 
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1999). Some common patterns, such as the "two degr ees separation" ([Williams et al. 



20021 ) and skewed or power law degree distributions (jDunne et all 120021 ). are 
found on binary food webs. However, these studies never considered weights, 
which may hide important informatio n of energy flux transferred between 
different species rtZhang and Guol . |2010| ). 

Weights on edges and nodes stand for energy-ma tter flux between tw o 



compartment s and biomass on each unit r e specti vely flHigashi et all Il993l ) 



Interestingly, iLindemanl (119421 ) and lOdum (ll988lV s sem inal works on food 
webs are all based on weighted networks. iPattenl (119851 ) et al. further de- 
veloped a systematic method called environ analysis to uncover some hid- 
den information in energy flows on networks. Besides the common phenom- 
ena found in earlier literatures inclu ding the hierarchical trophic str ucture 



and pyramid of biomass distribution (jOduml . Il983l ; lOdum et all 12004). some 



quan titative and ubiquitous patterns such as dom inant indirect effects (IHigashi et al 



19901 ). pathway prolife r ation (IPattenl . Il985c iBorrett et all 120071 ) and net 



1993), network amplification dPatten et all Il990f), network homoge nization (jPatten et al. 



work synergism (IPattenl. Il992f) are found by environ analysis in various food 
webs dFath and Patter] . Il999h . 

Metabolic theory is one of the greatest progresses in ecology in recent 
years which can be incorporated into food web studies (jBrownl . 120041 ) . Sev- 
eral universal pa tterns or laws re l ated to body size are d i scove red in the 



last two decades (IBrown and Westl . |2000| ; IWest and Brownl . 120051 ). For ex- 



ample, the three quarters power law relationship between metabolism and 
body mass (Kl e iber's law) is one of the most fundamental laws in metabolic 



theory (iKleibeil Il932h. Some ecologists also tried to link these patterns to 



trophic structures dNee et allll991 ;lLoeuille and Loreaul . 12006c iJennings and Mackinsonl . 



2003; iDamuthl. Il98ll; [Allen et all 120021). includin g the energetic equivalence 



rule ( Allen et al. , 2002 ; Loeuille and Loreaul . 2006 ); the trivariat e relat ionship 
among body mass, trophic level and energy flows (j Cohen et all 120031 ) . Nev- 
ertheless, some simple but important relations, say, biomass v.s. throughflow 
of each species, throughflow v.s. input flow and output flow are seldom ad- 
dressed by these previous studies. 

In parallel with these studies in ecology, the complex network commu- 
nity a lso started to pay attention to weight information of networks in recent 



years(lBarrat et all l2004j ; lAlmaas et all 120041 ; ISerrano et all 120091 ; iTumminello et al. 



20051 ). By incorporating the statistical mechanics method and random graph 



theory, weighted network analysis also revealed a series of universal patterns 



in various weighted networks, e.g. air traffic networks ( IBarrat et all 12004 ; 
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Guimera et all 120051). metabol ism networks (I Almaas et al.l . 12004 ). world trade 



web(Bhattacharya et all 120071 ) and stock-sharing networks of companies tlVitali et al. 
201ll ). etc. The new found common patterns include: long tailed distribu- 



tion of node intensity (total weig hts of each node), th e power law relation- 
ship between degree and intensity (IBarrat et all 120041) . the so called g ravity 



law (I Anderson! . 1201 ll ; lErlander and Stewartl . Il990l ; iKrings et al.l . 120091 ). and 
so forth. Energy flow networks in ecosystem no doubt are also weighted 
networks though the weight here h as the special meaning , i.e. the energy 
flux transferred by different species tjZhang and Guol . l2010l ). Therefore, it is 
reasonable to conjecture that the patterns found in other weighted networks 
should be also suitable for the weighted food webs. This paper tries to apply 
the approaches developed by complex weighted network studies to the energy 
flow networks in ecology. 

This paper is organized as follows: in Section [2j we introduce some ba- 
sic variables including the weights of edges and nodes. And also, the so called 



DGB D (Discrete version of a Generalized Beta Distribution) (IMartnez-Mekler et al. 



20091 ) curve which can fit the weights distributions better than the traditional 
curves is introduced. After that, the results of biomass and energy flow dis- 
tributions and several universal relationships including the allometric law at 
the population level, the dissipation law and the gravity law on 19 weighted 
food webs are shown in Section [3j After that, several interesting problem 
around these common patterns are discussed. We found an interesting neg- 
ative linear relationship between exponents of fitted DGBD distributions in 
energy flow, biomass and degree distributions (see Section 14. ip . And the 
mathematical relations among the scaling exponents are derived. Finally, 
the connection with the abundance-body mass relationship is discussed in 
Section HI 



2. Materials and Methods 

2.1. Data source 

We have investigated 19 food webs in different ecological environments. 
The food web information includes node (species or non-living compartment), 
node weight (biomass of a node), edge (energy flow relationship but not feed- 
ing relationship), and edge weight (the amount of energy flow from node i to 
j). The energy flow between two nodes was measured as the unit volume flow 
of the carbon element into or out of the node (the unit is gC/m 2 /year)). The 
biomass stands for the total mass of living biological organisms of a species in 
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Table 1: Empirical food webs and their topological properties 
(N stands for the number of vertices of the network and E is the number of 
edges. The webs are sorted by E.) 



Food web 


Abbre. 


N 


E 


Crystal River Creek (Delta Temp) 


CrystalD 


23 


60 


Crystal River Creek (Control) 


CrystalC 


23 


81 


Chesapeake Bay Mesohaline Net 


Chesapeake 


38 


122 


Lower Chesapeake Bay in Summer 


ChesLower 


36 


115 


Middle Chesapeake Bay in Summer 


ChesMiddle 


36 


149 


Upper Chesapeake Bay in Summer 


ChesUpper 


36 


158 


Narragansett Bay 


Narragan 


34 


158 


Lake Michigan 


Michigan 


38 


172 


St. Marks River (Florida) 


StMarks 


53 


270 


Mondego Estuary - Zostrea site 


Mondego 


45 


348 


Cypress, Wet Season 


CypWet 


70 


545 


Cypress, Dry Season 


CypDry 


70 


554 


Everglades Graminoids, Dry Season 


GramDry 


68 


793 


Everglades Graminoids, Wet Season 


Gram Wet 


68 


793 


Mangrove Estuary, Dry Season 


MangDry 


96 


1339 


Mangrove Estuary, Wet Season 


MangWet 


96 


1440 


Florida Bay, Wet Season 


BayWet 


127 


1938 


Florida Bay, Dry Season 


Bay Dry 


127 


1969 


Florida Bay 


Florida 


127 


1938 



a certain period of time and per unit volume. Customarily it was also mea 
sured by carbon content (the unit is gC/m 2 ) (IBaird and Ulanowiczj. Il989l ) 



These food webs' information is obtained from the online database 1 , which 



is ba s ed on the published papersdBaird et all Il998l; IBaird and Ulanowicz , 



19891: llJlanowicl Il986l : lAlmunia et all Il999l : iMonaco and Ulanowiczj . 119971 ; 



Hagy . 20021 ) . In Table [H we list the name and the number of nodes N and 



edges E of each web. 

2.2. Basic variables 

Our work is based on the flux matrix of a weighted food web. An eco- 
logical energy flow network is a weighted directed graph that represents re- 



http://vlado.fmf.uni-lj.si/pub/networks/data/bio/foodweb/foodweb.htm 
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Figure 1 : An Illustration of variables 



lationships of ecological energy transfer between species. This graph can be 
represented by a flux matrix: 



\N+2)x(N+2) 



{fij}(N+2)x(N+2), Vz,j G [0,JV+1] 



(1) 



where f\j is the energy flow from species i to j. Two special nodes rep- 
resenting environment (node and node N + 1) are added to the web. 
Node denotes the source of energy flow, whereas node N + 1 represents 
the sink. The dissipative and exported energy flow to the node N + 1. 
Therefore, there are totaly (N + 2) x (N + 2) entries. The flux matrix 
/r;v+?w(Ar+2) can be read from the original weighted food webs(IBaird et al 



1 99a iBaird arid Ulanowiczl . Il989t lUlanowiczl . Il986l : lAlmunia et "all . Il999l ; | Monaco and Ulanowiczl . 
1997tlHagyL l2002h . 



We can calculate the total through flow of any given node % according 
to the flux matrix F(tv+2)x(v+2) (see figured!) . This value is also c alled node 
strength in complex weighted network studies (jAlmaas et all 12004 ) . Because 
the network is always balanced, we need only to calculate the efflux of each 
node as 71- , 



N+l 



N 



(2) 



j=0 



In addition, we define another variable to represent the weight of a node 
Xi, indicating the biomass of i. This information is also available from the 
original weighted food webs. 
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Di is the dissipated flow to the environment (node N + 1) from species 
v. fi t N+i ( see figured]). The dissipation flow includes the flows of output and 
respiration. 

kf 1 and k° ut are in-degree and out-degree, i.e., the number of inward edges 
(excluding edges from the source) and outward edges (excluding the edges to 
the sink) of i respectively. In the example network in figure [U kf 1 = 3 and 

k out = L 

2.3. Distributions and DGBD curves 

We will study the distributions of Tj, Xj, k\ n and k° ut in any empirical 
food web. Instead of g iving the empirical density or distribution function 



(IZhang and Guol . 120101 ). we use the rank-ordered curve to show the distri- 
butions of these variables. For example, if we have a small food web with 5 
species, and their biomass values are {100, 19, 200, 5, 1} gC/m 2 . A sequence 
of biomass values in a decreasing order can be obtained: {200, 100, 19, 5, 1}. 
Then we plot this sequence on a coordinate with the horizontal axis as the 
rank value, namely {1, 2, 3, 4, 5} and the vertical axis as the biomass values. 
So the final curve on this coordinate is the rank-ordered curve. The main 
advantage of adopting this curve is its simplicity for the calculati on, as well 



as it contains the same information as the distribution function (INewman 



2005|). 

Then, we use the DGBD (Discrete version of a Generalized Beta Distri- 
bution) function to fit the rank-ordered curve: 

Y{ n ) = A {N + l ~ n) \ a,b>0 (3) 

i 

where, Y(ri) is the value of the concerned variable (Tj or XA of the node 
i, and r« is the decreasing order of i ranked by Y(ri) values. N is the total 
number of nodes in the food web. A, a, b > are parameters to be esti- 
mated. A stands for the magnitude of flow or biomass in this food web 
which is dependent on the measurement units, a and b are exponents of 
power laws in the tail and head part of the curves respectively. If we set 
a = 0, then formula (j4]) becomes K(rj) = A/r%, which is the famous Zipf's 



law (INewmanl . 120051 ). However, the classic Zipf's law is not always the best 



choice for fitting empirical data due to the large deviat ions in the tail part 



of the rank-ordered curves (IMartnez-Mekler et all |2009i ). This disadvantage 
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can be mended by introducing a new exponent a in the DGBD fitting. Pre- 
vious study shows that form ula (jlj) can fit lots of empirical data very well 



flMartnez-Mekler et all 120091 ) 



Notice that, if we set a = in equation HI then the rank-ordered curve of 
a variable Y follows Zipf's law with the exponent b. Therefore, the density 
function of the random variable Y is a power law with the exponent 1/6+1 
due to the one-one corresp ondence between the rank-ordered function and 
probability density function ( Newman . 20051 ). So, DGBD curve is capable to 
not only fit those data with power law tails (note that the tail part of the 
density function is just the head part of the rank-ordered curve) but also the 
data with an obvious deviation from power laws by tuning the extra exponent 
a (figure [2]). 

As shown in figured different combinations of a and b correspond different 
shapes of the distribution curve. Exponent b is the slope of the head part of 
the rank-ordered curve in figure [2] (b) as well as the tail part of the probability 
density curve in figure [2] (c). Therefore, we say b indicates the heterogeneity 
of Y distributing for the large species. On the other hand, a indicates the 
unevenness of the tail part of the rank-ordered curve (figure [2] (a) and (b)) 
as well as the head part of the probability density curve (figure [2J (c)). As 
a increases, the tail of the rank-ordered curve drops very fast so that the Y 
sharing heterogeneity in the small species is very large. 

We mainly adopt OLS (Ordinary Linear Square regression) method to fit 
DGBD curves and other power law relationships. We can take logarithmic 
on equation [4] to derive: 



log(y(r<)) = \og(A) + a\og(N + 1 - r<) - 61og(n). (4) 

Thus, log(F(rj)) depends on two variables log(iV + 1 — ri) and log(rj) 
so that the bivariate linear regression method can be used to derive the 
coefficients log(A),a and b. 



3. Results 

3.1. The DGBD distributions 

We obtain the distribution of Tj, JQ, k™ and k° ut for each of the 19 em- 
pirical food webs, and fit them by DGBD curves. One selected food web as 
an example is plotted in figure [3j 

From figure [31 we know the curves can be divided into three parts by two 
inflexion points, and each part obeys independent logarithmic decreasing 
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Figure 2: Different shapes of rank-ordered curves with different a and b under the same 
A = 1 and TV = 100. (a) Rank-ordered curves on a linear-log plot; (b) The same rank- 
ordered curves on a log-log plot; (3) The probability density curves of the same four groups 
of data on a log-log plot 
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behavior. Obviously, the head and tail parts of the curve have much steeper 
slopes than the middle part. The local slope of the curves indicates the 
heterogeneities of the throughflow distribution. In other words, the larger 
the absolute value of slope is, the higher the degree of heterogeneity of a 
vertex correspondingly is. Therefore, we can conclude that the nodes in the 
heads and tails are more heterogeneous than the ones in the middle. 

In figure [31 we distinguish nodes by their trophic levels. Green circles, 
blue squares, and red triangles represent the first, second and third trophic 
level species respectively. It is observed that nodes of first trophic level locate 
both at the head and tail parts, while most of the second level species locate 
at the middle part, and most third trophic level species are in the middle 
or tail parts of the curve. This distribution pattern of species on different 
trophic levels is similar for all large food webs (the food webs below Michigan 
web in tabled]). We may conclude that the distribution of throughflow on the 
second trophic level is much more even than the first and third trophic levels. 
The inset figure of figure [3] plots the rank-ordered distributions of in-degree 
and out-degree of all edges. These two curves have the similar shape with 
the energy flow and biomass distributions. 

In Table [21 we list all the fitted parameters and R 2 s of DGBD for the 19 
webs. Only the exponents of a and b are shown in the table because A repre- 
senting the magnitude of throughflow or biomass is relatively unimportant. 
By comparing different rows, we know that the food webs with more edges 
can be better fitted by DGBD curves because their R 2 s are larger. Notice that 
there are three food webs (ChesLower, ChesMiddle and CheslJp) having ex- 
ponents = and one web (CrystalD) having exponents ax t , a k in, a k in = 0. 
That means these webs are anomalous and reach the extremal cases of DGBD 
curves. However the R 2l s of the the fitting curves are also very high. 

The last four columns show the DGBD fittings of in and out degrees. 
An obvious trend is the exponent a ki increases, however b ki decreases with 
the size of the food web. That indicates the number of connected species 
distributes on different nodes more evenly in small food webs and the average 
degrees increase with the size of the web. 

3.2. Allometric scaling laws at the population level 

According to the similarity between the distribution curves of Tj and Xi in 
figure [31 we guess there may exist a connection between variables Xi and Tj. 
The log-log plot of these two variables shows that the relationship between 
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Figure 3: Rank-ordered distributions and DGBD fittings of Tj,X;, k° ut and k\ n for Mang- 
wet food web. The indications of 1st, 2nd and 3th level species in the legend represent 
the first, second and third trophic level species in the food web. 
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Table 2: DGBD fitting parameters for different variables 
(The webs are sorted by their number of edges.) 



TTT 1 

Web 


a Tl 


7 

t>Ti 


D2 

itrp 




7 


D2 

K Xi 




7 

Ohout 




7 

Ufein 


CrystalD 


r\ at* 

0.45 


3.96 


0.96 


0.00 


A. 21 


0.97 


0.00 


0.69 


0.00 


0.98 


CrystalC 


0.76 


3.41 


0.97 


0.24 


3.89 


0.95 


0.20 


0.56 


0.35 


0.71 


Chesapeake 


2.93 


1.54 


0.94 


1.70 


2.17 


0.99 


0.23 


0.50 


0.19 


0.73 


ChesLower 


10.02 


0.00 


0.80 


5.53 


0.75 


0.85 


0.68 


0.17 


0.15 


0.88 


ChesMiddle 


6.88 


0.00 


0.83 


3.42 


0.77 


0.87 


0.55 


0.25 


0.40 


0.63 


ChesUpper 


4.77 


0.00 


0.72 


3.42 


0.86 


0.88 


0.61 


0.16 


0.37 


0.70 


Narragan 


0.80 


2.54 


0.93 


0.16 


1.88 


0.93 


0.21 


0.47 


0.41 


0.56 


Michigan 


6.59 


0.81 


0.95 


4.80 


0.18 


0.89 


0.51 


0.47 


0.45 


0.45 


StMarks 


1.25 


1.33 


0.99 


0.94 


2.14 


0.98 


0.59 


0.38 


0.54 


0.40 


Mondego 


3.36 


1.70 


0.97 


1.41 


2.43 


0.99 


0.57 


0.26 


0.48 


0.94 


Cypwet 


2.74 


2.63 


0.97 


3.03 


4.23 


0.98 


0.41 


0.56 


0.62 


0.66 


Cypdry 


2.34 


2.36 


0.95 


2.51 


3.99 


0.98 


0.43 


0.52 


0.63 


0.66 


Gramdry 


2.08 


3.18 


0.97 


2.75 


2.75 


0.97 


0.34 


0.43 


0.67 


0.54 


Gramwet 


3.03 


2.94 


0.97 


3.74 


2.40 


0.97 


0.34 


0.43 


0.67 


0.54 


Mangdry 


1.51 


3.22 


0.97 


1.41 


4.08 


0.98 


0.66 


0.32 


0.79 


0.42 


Mangwet 


1.67 


3.32 


0.98 


1.63 


4.18 


0.99 


0.66 


0.32 


0.79 


0.43 


Baywet 


1.65 


2.80 


0.97 


2.26 


2.40 


0.99 


0.74 


0.40 


0.97 


0.25 


Baydry 


1.62 


2.63 


0.98 


2.24 


2.28 


0.98 


0.75 


0.39 


0.97 


0.25 


Florida 


1.65 


2.80 


0.97 


2.26 


2.40 


0.99 


0.74 


0.40 


0.97 


0.25 
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Table 3: Fitting exponents and goodness of power law relationships 



( The webs are sorted by their number of edges.) 



Web 


a 




R 


2 

alio 


c 









d2 

^diss 


7 




d2 
-^qra 


7i 


72 


R 


2 

hi 


V 


CrystalD 





98 


0. 


90 





67 





96 


1 


00 





63 





70 





57 





75 


0. 


74 


11.5 


Crystal C 





92 


0. 


88 





68 





96 





99 





53 





65 





50 





57 


0. 


65 


2.125 


Chesapeake 


1 


03 


0. 


75 





50 





99 





98 





68 





84 





62 





77 


0. 


85 


-9.33 


ChesLower 


1 


75 


0. 


91 





49 





95 





99 





70 





75 





61 





84 


0. 


76 


-1.33 


ChesMiddle 


1 


58 


0. 


82 





47 





88 





85 





67 





77 





60 





78 


0. 


78 


-1.43 


ChesUpper 


1 


21 


0. 


64 





58 





95 





99 





64 





64 





62 





67 


0. 


64 


-2.19 


Narragan 


1 


17 


0. 


55 


2 


05 





81 





94 





54 





81 





49 





60 


0. 


81 


-2.47 


Michigan 


1 


20 


0. 


94 





67 





99 


1 


00 





62 





86 





57 





72 


0. 


87 


-2.25 


StMarks 





69 


0. 


70 





43 





99 





95 





68 





74 





76 





56 


0. 


75 


-0.19 


Mondego 


1 


24 


0. 


87 





53 





98 


1 


00 





79 





85 





83 





70 


0. 


86 


-2.04 


Cypwet 





66 


0. 


78 





46 





97 





99 





70 





84 





85 





55 


0. 


87 


-0.26 


Cypdry 





63 


0. 


76 





41 





96 





95 





68 





81 





81 





57 


0. 


83 


-0.32 


Gramdry 





90 


0. 


90 





58 





97 


1 


00 





66 





76 





61 





73 


0. 


77 


1.5 


Gramwet 





93 


0. 


92 





59 





98 


1 


00 





71 





81 





66 





79 


0. 


81 


2.57 


Mangdry 





77 


0. 


82 





45 





98 





98 





58 





77 





60 





56 


0. 


77 


0.09 


Mangwet 





78 


0. 


83 





44 





98 





98 





59 





77 





60 





57 


0. 


77 


0.13 


Baywet 





87 


0. 


81 





33 





92 





95 





62 





79 





67 





54 


0. 


80 


0.92 


Bay dry 





85 


0. 


81 





32 





91 





95 





61 





78 





68 





52 


0. 


78 


0.67 


Florida 





87 


0. 


81 





33 





92 





95 





62 





79 





67 





54 


0. 


80 


0.92 



Xi and Tj is actually a power law: 



Ti oc XT 



(5) 



where exponent a is a parameter to be estimated for each food web. As 
shown in upper row in figure HI the sample points aggregate around their 
fitted lines very well. This relationship is ubiquitous for all 19 food webs as 
shown in Table [3j We use the ordinary linear regression method to find the 
best fitting lines (Table |3]). And most _R 2 's are larger than 0.8. The i? 2, s 
and exponents decrease with the scale of the network because the statistical 
significance increases with the number of samples. 

This specific scaling relationship reminds us the famous allome t ric scal- 



i ng law (jBrownl . 12004c IWest and Brownl . 120051 ; iBanavar et al.l . Il999l ) . iKleiber 



(119321 ) found that the metabolism and body size of an organism usually fol- 
lows a ubiquitous power law relationship with an exponent around 3/4. If 
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Figure 4: Power law relationships of Tj v.s. Xi and Di v.s. JQ. The original data points 
as well as the OTS fittings are shown on the log-log coordinate 
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we treat the whole population of a species as an integrated organism, then 
Tj is its metabolism and Xi is its body mass. Therefore, equation can be 
viewed as the allometric scaling law at the population level. Nevertheless, 
unlike the universal Kleiber's law for species, the allometric scaling expo- 
nents of population on food webs are not universal but fluctuate in between 
[0.63,1.75]. 

3.3. Dissipation law 



As pointed by the earlier ecological studies ( lOduml . Il983l ; iLindemanl . Il942[ ) . 
a large fraction of energy flows dissipates to the environment in the whole 
ecosystem. The dissipated energy flow can be captured by the variable Di 
which is also available from the original data. Empirically, scales with 
throughflow Tj in the following way: 

A = elf . (6) 

Equation is called dissipation law in this paper, where c and /3 are param- 
eters to be estimated. We observe that the estimated exponents (3 are all 
slightly smaller than 1 (see the 2nd row of figure H] and the 4th 6th columns 
of table |3]), therefore the dissipation rate (dissipation per throughflow) de- 
creases with the throughflow of the species slightly. If (3 = 1, c is the average 
energy dissipation rate of the whole food web. Since the empirical exponents 
in table |3] are approaching to 1, the coefficient c's are almost the dissipation 
rate of the specific food web. From table [31 we can read c's are fractional 
numbers that are smaller than one except the food web Narragan whose 
exponent (3 deviates 1 significantly. 

3.4- Gravity Law 

Although we have studied the ways of energy flows correlate with the 
biomass and dissipation, we still don't know how the energy flows distribute 
among different pairs of species. Actually, the large throughflow nodes can 
exchange large energy flows. This effect is reflected by the so called gravity 
law, namely, the energy flow between i and j scales with the product of the 
total throughflows of i and j, i.e., 

Ju oc (Til))'. (7) 

This scaling relation is known as the gravity law in other complex systems. 
Researchers found the flows, say traffic flow or trade flow between cities 
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(a) (b) 
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Figure 5: Univariate and bivariate gravity law of Mangwet food web. (a) The uni- 
variate scaling relationship fa oc (TiT)) - 587 ; (b) The bivariate scaling relationship 
fij oc t°- 6 Tj- 57 . The original data points as well as the OLS fittings are shown on the 
log-log coordinate 

or countries scales with {mxm%)' 1 / <F ^ where m« is the size of the system 
(population of a city or GDP of a country) and d is the distance between the 
two systems, 7 and r are fitting exponents. In our case, the throughflow of 
each node Tj is treated as the size of a node comparable to the population 
of a city. However, the distance d has no correspondence because spatial 
information is not included in our food webs. 

Figure |5^a) shows this phenomenon and the parameters are listed in the 
last two columns of table EJ This pattern is not as significant as the previous 
two patterns because the R 2, s are always smaller than 0.85. From figure [5] 
(a), we could observe that there are several straight bands in the clusters of 
data points which indicates that the energy flow between two nodes may 
be predicted by other variables rather than the product of Tj and Tj. 

As the studies in gravity laws, we suggest that the following bivariate 
scaling relation holds, 
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fu Ty. (8) 

Where, 71 and 72 are two estimated parameters. Figure |5] (b) shows 
this bivariate gravity law. All the estimated exponents 71, 72 and the cor- 
responding R 2 are shown in the last 2nd, 3rd and 4th columns in table El 
The bivariate gravity law equation |8] can fit the original data better than 
equation [7] since the corresponding -R^'s are slight higher than i?L a 's. An- 
other interesting phenomenon observed from table [3] is that the exponent 72 
is larger than 71 for small food webs but smaller than it for large food webs 
except Gramdry and Gramwet. We can conclude that the large energy flows 
prefer to link nodes with large throughflow in all food webs. The dependence 
of energy flow on the donator (prey) and receptor (predator) is asymmetric. 
As the energy flow along each edge increases, the energy throughflows of 
receptors increases faster than donators in small food webs. The speed of 
increasing of throughflows for donator is higher than receptors in large food 
webs. 

All of these observed patterns of energy flows exhibit statistical signifi- 
cance and universality for all 19 empirical food webs. 



4. Discussions 

4-1. Relationships of exponents 

In section 13.11 and table [2J we have obtained a set of fitting exponents of 
DGBD curves. These exponents also have some patterns. It is observed from 
table [2] that exponents a and b have a negative correlation. The relationship 
is very clear once we plot the pairs of (a, b) in one coordinate (see figure [6]). 
We separate the a, b exponents of Tj, Xj and the ones of k° ut , k\ n because the 
former has a wider range. However all these pairs of (a, b) show the nearly 
linear relationship with similar negative slopes (the mean slope is —0.5) and 
statistical significance (the R 2, s of these relationships are all larger than 0.5 
except kf 1 ). 

We suppose this pattern reflects a kind of particular regularity o f food 



webs since the similar phenomenon is never reported in previous studies tlMartnez-Mekler et al. 



20091 ). The negative correlation between exponents a and b implies a kind 
of complementarity between the heterogeneities of energy flow or biomass 
resources distributions at the head and tail part species on rank curves. As 
the heterogeneity of the small species increases, the unevenness of energy 
distribution in large species decreases. 
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Figure 6: The linear relationships of DGBD fitted exponents a and b in four distributions 
(T. l ,X l ,k° ut ,kl n ). The R 2, s of the best fitting lines are R\ = 0.51, R\. = 0.72,i^ ?ttt = 

0.67, R\ m = 0.35 
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(the 45 degree diagonal straight line). 
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Besides this pattern, other two interesting relationships between the ex- 
ponents of distributions and power law relationships can be derived . We can 
write down the DGBD distributions of Tj and Xi in the following forms: 



r. 



and 



x , = A x <» + 1 -«r\ (10) 



for each node i. As we have shown in section 1372] Tj and Xi follow a power law 
relationship: equation (JSJ). If we insert equation (J5J) and ffTOj) into equation 
(E]), we can easily derive the following equation: 

^ bx, ' ^ ' 

i i 

This should be satisfied for any r,. So comparing the coefficients of the terms 
(N + 1 — Tj) and r i; we can derive the following relationships. 

a T% = a x ,a, (12) 

and 

b Ti = b Xi a. (13) 

If these two relations hold for all food webs, the pairs of (a^aaxj or 
(pTiicAxi) in the empirical food webs should form a straight line with a 
45 degree slope. 

From Figured the pairs of (ay^aaxj or {br^abx^) concentrate around 
the predicted relationship with small deviations. That means the predicted 
relationships (equations [12] and [T3]) are almost correct for the empirical 
food webs. However, a systemic deviation from the predicted line exists 
for (ciT^ctciXi) since all the data points are lower than the theoretical line. 
Therefore, we may underestimate the values of ax, or overestimate a^- The 
errors can not come from the estimation of a because the similar systemic 
deviation of is not observed. However the reasons for these errors are 
still mysteries for us. 
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4-2. Other possible patterns 

Besides the common patterns shown in the previous texts, we have also 
investigated other possible patterns exhaustively. However they are either 
unclear or trivial. 

For example, there is a power law relationship between node degree and 
strength in other w eighted complex networks as shown in previous studies 



( iBarrat et all 12004 ). Although a positive correlation between the degree 



and strength can be observed in our food webs, this power law relation is not 
significant (R 2 = 0.22). 

The energy flow distributions for each node also obey DGBD curves, but 
this content is abandoned in the main text because it can not provide us 
more insights. 

Another trivial scaling relationship is between the energy flow and the 
product of biomass of the two species XjXj because it is an obvious result 
from the scaling relationship in equation [7] and equation 

Finally, one may guess that a power law relationship between ki n and k out 
must exist because of the similarity of the distribution curves in the inset of 
figure |2j However, their relations are not significant once we draw the pairs 
of ki n , k ou t in one coordinate. Therefore, all the patterns we have selected are 
significant and nontrivial. 

4-3. Abundance and Body Size 

Body size is treated as a very fundamental observable in ecology because 
it determ ines other important variables of organisms including the t rophic 



struct u re (ICohen et all l2003l ; lBrownl . 12004 ; iBrown and Gilloolyl . 120031 ) . ICohen et al 
(l2003h . lBrown and Gillooivl fl2003l ) discussed the scaling relationship between 
abundance and body mass in food web. We may derive a relationship between 
abundance and body size from the scaling relation between throughflow and 
biomass of each node. Suppose the numeric abundance of species i is N i: 
the average body mass of % is Mj, and its average energy metabolism is Fj. 
Then, according to the definitions, 



and 



Xi = NiM i} 



F = NiFi. 



(14) 
(15) 



The Kleiber's law links Mj and Fi as follows for any species i, 
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3/4 



Ft oc M\ 

So, insert these relations into equation EJ we have, 

3/4-g 

Ni oc Mi a ~" . 



(16) 



(17) 



Therefore, the scaling exponent between abundance and body size r\ is: 

3/4 -a 



a — 1 



This exponent can be estimated according to the exponent a for any em- 
pirical food webs. The rj values are derived in the last column of Table |3j 
We can see that the exponent a may be either positive or negative which 
means the abundance may increase or decrease with body mass. This con- 
clusion contradict s with our ob s ervat i ons and the exponent devia tes from the 
previous studies (j Cohen et al.l . 120031 ; iBrown and Gilloolyl . 120031 ). We guess 
the problem is the nodes in our webs do not always stand for living species 
but other non-living compartments, therefore equation [T7] may not hold for 
all nodes in the whole network. As a result, the abundance-body size scal- 
ing exponent can not be determined by a solely. More discussions on linking 
trophic structure, energetics and metabolic theory are deserved for the future 
studies. 



5. Concluding Remarks 

The weighted food webs have several common patterns that have never 
been found in previous binary food web studies. First, the energy flow, 
biomass and degree resources distribute on different species far from evenness. 
This heterogeneity can be characterized by a common rank-ordered curve 
called DGBD. Second, there are a set of scaling relationships in weighted 
food webs. The power law relationship between Tj and X{ can be regarded 
as a counterpart of Kleiber's law in population level. The scaling relation be- 
tween the dissipation and throughflow characterize the dissipative nature of 
energy flows in ecosystems. And another interesting common pattern is the 
so called "gravity law" which is also discovered in other complex systems. 
Finally, we find two interesting regularities in the fitted exponents includ- 
ing the negative correlation between a and b and the predicted relationship 
between the distribution exponents and the power law relation exponent. 
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This paper only exhibits these common patterns in the empirical weighted 
food webs, however the underlying mechanisms that can reproduce these 
patterns are left for the future works. 
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